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Abstract 

The corrections to neutrino mixing parameters in the presence of matter of constant density are 
calculated systematically as series expansions in terms of the mass hierarchy Arn^i/ i^m'^^. The 
parameter mapping obtained is then used to find simple, but nevertheless accurate formulas for os- 
cillation probabibilities in matter including CP-effects. Expressions with one to one correspondence 
to the vacuum case are derived, which are valid for neutrino energies above the solar resonance 
energy. Two applications are given to show that these results are a useful and powerful tool for 
analytical studies of neutrino beams passing through the Earth mantle or core: First, the "dis- 
entanglement problem" of matter and CP-effects in the CP-asymmetry is discussed and second, 
estimations of the statistical sensitivity to the CP-terms of the oscillation probabilities in neutrino 
factory experiments are presented. 



1 Introduction 



With the development of long baseline neutrino beams passing through the mantle of the Earth, 
three flavor neutrino oscillation with a constant matter profile is presently drawing attention. Some 
effort has been spared on the exact solution of the connected cubic eigenvalue problem [|^. However, 
the obtained solutions are huge and were up to now only used in computer based calculations. 
Also approximative solutions for oscillation probabilities and mixing angles have been proposed 
for several parameter regions |2|, which are interesting and useful. The intention of this work 
is to first derive analytic approximations for the mixing parameters in matteij^ according to the 
standard parameterization, which then allows to compute all desired quantities like probabilities 
or amplitudes from the known expressions in vacuum by substitution. The parameters in matter 
are calculated in a series expansion in the small mass hierarchy parameter a := Arn^i/ Am'^^. The 
obtained results are discussed and then applied to the appearance channel probability P(i^e t^h)- 
A simple solution, which is easy to use, but nevertheless accurate over a wide parameter range is 
obtained. No new notation is introduced besides abbreviations known from two neutrino oscillation 
in matter. Furthermore, the result shows at first sight the convergence to the vacuum case at 
small baselines and thus is directly connected to the terms in vacuum. The approximate solutions 
obtained with this method are a powerful tool for further analytical studies. To demonstrate this, 
two applications are given. First the derived expressions are exploited to compute the frequently 
used quantity called the CP-asymmetry A^^ , which has considerable importance in CP-violation 
studies. The problem is that matter effects cause contributions to the CP-asymmetry, which 
cannot easily be distinguished from intrinsic CP-effects. Here, expressions for A^^ in matter are 
given for high neutrino energies (more precise: low L/Ey). The result is then used to investigate 
what can be learned from the energy dependence of A'^^ . The second application given estimates 
the statistical sensitivity to the CP-terms of the oscillation probabilities in neutrino factory long 
baseline experiments. Plots are presented, which show the magnitude of CP-effects at different 
baselines and beam energies. Contrarily to what presently can be found in the literature, the here 
obtained results indicate strongly that, in general, the low energy option is not the best solution 
to measure effects from the CP-phase 5. The reason for this discrepancy is discussed. 

Throughout this work two assumptions will be made: First, that the mass hierarchy parameter 
a := Amfg/Am^^, which is used as expansion parameter, is small. Consider for example an 
atmospheric /\m? of 3.2 • 10~^ eV^ Q . For solar mass differences of LMA-scale0 Q between 
10-5 eV^ and 10~'^eV^ a varies between 0.0031 and 0.031. Second, it will be assumed that 
the mixing angle ^13 is small as indicated by reactor, solar, and atmospheric experiments. The 
strongest bound is given by the CHOOZ experiment [|| with sin^ 2^13 < 0.1. The smallness 
of this parameter will be used to classify terms, which appear in the expressions for oscillation 
probabilities. The mixing angles 9i2 and ^23 should be chosen from the interval [0,7r/2]. 



^ Oscillation in matter can be described by a mapping of the six basic parameters 6^12, ^23, Am|i, Am|i, 
and 5 similar to the well-known two neutrino oscillation formulas in matter. 

^The abbreviation "LMA" stands for Large Mixing Angle MSW-solution to the solar neutrino problem. The 
MSW-solution assumes resonance enhanced oscillation of neutrinos passing the core of the sun. 
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2 Three neutrino oscillation in vacuum 

In vacuum, the neutrino oscillation probabilities are given by the well-known formulas 

i>j i>j 

with the abbreviations J^™ := UiiU*-U^^Umj and Ajj := AmjjL/{AE). Here, [/ is the mixing 
matrix of the neutrino sector in standard parameterization form: 

/ C12C13 C13S12 e"*'^si3 \ 

U = -S12C23 - e"^Ci2Si3S23 C12C23 - e*''si2Sl3S23 C13S23 • (2) 

\ -e*'^Ci2Si3C23 + S12S23 -e*''si2Sl3C23 " C12S23 C13C23 / 

Since in this work, the hierarchy |Am|;^| <C |Am|^| between the two mass squared differences 
is exploited, from now on all mass squared differences will always be related to the atmospheric 
squared mass difference: Am|-^ =: A, Am^i = aA, Amig = (1 ~ ct)^i and A = AL/(4£^). 
Series expansion up to order gives the following important terms in the oscillation probability 



Po = sin^ 6*23 sin^ 2(9i3 sin^ A (3a) 

-fsin (5 = a sin 6 cos ^13 sin 2^12 sin 2^13 sin 2^23 sin^ A (3b) 

-fcos 5 = a cos 6 cos 013 sin 2^12 sin 2^13 sin 2023 cos A sin^ A (3c) 

P3 = cos^ (923 sin^ 2(9i2 sin^ A (3d) 



Expanding the oscillatory terms in a means linearization of the oscillation over the solar mass 
squared difference. This gives valid results only for a A < 1. With todays knowledge about 
neutrino masses this does not cause crucial errors for neutrino energies above 1 GeV at baselines 
below approximately 10000 km. The two terms PsmS and -Pcos^i containing the CP-phase 6, are 
both of order a and hence suppressed by the mass hierarchy. This reflects the fact that CP-effects 
vanish when the mass hierarchy becomes large. Besides the factor sin^ 623, the term Pq is similar to 
the two neutrino oscillation probability which in matter is expected to show the resonant behavior 
called MSW-effect |^]. The term P3 is the only term of order a^, which is not suppressed by the 
small mixing angle ^13. Hence, it is important to take this term into account when ^13 is small. 
If 013 is not too far away from the CHOOZ-bound, P3 can safely be neglected. All other terms of 
order are additionally suppressed by one or more powers of 0i3 and are not listed here. 



3 Mixing parameters in matter 



In matter, the effective Hamiltonian in flavor basis is given by 
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(4) 



Here U = ?723(023)f^i3(0i3) '^)t^i2(0i2) is the mixing matrix, which rotates from mass to flavor basis. 
The second term is generated by matter effects with A = 2VEy and V = \f2Gpne-, where Gi? is 
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the Fermi coupling constant and rif. is the electron density of the matter, which is crossed by the 
neutrino beam. 

The matter term is invariant under rotations in the 23-subspace. Separating diag(mf , mf , mf ) 
which, as global phase, does not contribute to the probability, and using the above defined para- 
meters, the Hamiltonian can be written in the form 
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(6) 



the relations 



ulUni0i3,S)Us 

ului2{ei2)us 

ut diag(o, b, c) Us 



f/l3(%,0) , 
Ul2{9l2) , 

diag(a, b, c) 



are valid. Inserting the identity matrix Us U^ at the appropriate places in eq. (^) gives 



n = ^U2sUs 



Ui3{ei3,o)u 
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Diagonalization of the real matrix Af by U := ^^23(^23) ^^13(^13) ^^12(^12) together with the part 
which was factored out gives the complete mixing matrix U' in matter: 



U' = ?723(^23) Us U2z{e2z) t/l3(%) t/l2(^12) ■ 



(9) 



Mixing angles in standard parameterization form 

The matrix U' must still be brought to the standard form. The matrix 



^723(^23)t/5t/23(^23) 



1 

c s 

-e*'^^* e^^C* 



(10) 



with 



C := cos ^23 cos 



^23 - sin 6*23 sin ^23 



S 



cos 6*23 sin 6*23 + e sin 6*23 cos 6*23 



(11a) 
(lib) 
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can be made real by the phase rotations /3 := — argC, 7 := argS*, and 5' := argC — argS* | 



10 o\ /100\/100\ 

e-'f U23{e23)UsU23{023 ) 1 = |C| \S\ ] . (12) 

-e(*^-T) / V e-'^' J V \C\ J 

This gives 

/lO \ /io 0\ 

U'=\0 e'f |C| |S| U/5'[/l3(^13)C/5Wl2(^12)C/5' . (13) 

V -e(*^-^) / V \C\ J 

The phase rotations on the left and on the right can be absorbed in the field vectors, yielding then 
U' in standard parameterization form: 

U' = U{e'23) UisiOvi, 6') [/i2(0"i2) . (14) 

This finally means, that the (standard) mixing angles and in matter are equal to ^13 and 
9i2 which are obtained from the matrix that diagonalizes M. The matter correction 823, however, 
mixes with the CP-phase 6: 

e'u = 9i3, (15a) 
0'i2 = Ou , (15b) 



2/ 22^ 22^ ^ / \ 

sin ^23 = cos 6*23 sin ^23 + sin ^23 cos ^23 + 2 cos 5 sin ^23 cos 6*23 sin 023 cos 6*23 , (15c) 

. , . sin 2023 f.^.. 

smo = smo^ — -— — . (15dj 



23 



Equation (|l^d) was first found by S. Toshev There, a different parameterization is used, which 
- for oscillations - is equivalent to the standard parameterization. It is important to note that 
the results given up to here are exact results for three neutrino oscillation in matter and do not 
presume that the mass hierarchy parameter is small. 



Calculation of the eigenvalues and eigenvectors 

Hereafter A will be used as abbreviation for ^. Diagonalization of the matrix M leads to the 
oscillation parameters in matter. Note that M does not include the parameters 023 and 5, which 
have been factored out. This will simplify the calculation of the eigenvalues and eigenvectors of 
M considerably: 



■S13 + ^ + acfgsfa asi2Ci2Ci3 S13C13 - asi3Ci3sf2 

M = I aSl2Cl2Cl3 Qcf2 -aSl2Cl2Sl3 | ■ (16) 

S13C13 - aSl3Cl3'Si2 -aSl2Cl2Sl3 + asf2S?3 



Using \C\ and \S\ in eq. ( |l2[ ) further restricts the parameter space for 623. Since 623 is assumed to be close to 
7r/4 and 623 in general is small, this problem is not relevant for the calculations presented here. 
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The invariants of the cubic eigenvalue problem are given by 

h = Tr(M) = Ai + A2 + A3 = 
= A+l+a , 

h = \ [TV(M) - Tr(M2)] = A1A2 + A1A3 + A2A3 = 

= A cos^ ^13 + a + a A (sin^ ^13 sin^ 612 + cos^ ^12) , 
73 = Det(M) = A1A2A3 = 
= cos^ ^13 cos^ ^12 . 
Solving this system of equations in a series expansion of a gives the eigenvalues 



Ai 
A2 
A3 



2^ ' 2C 

2^ ' 2C 



with 



C =\l{A-cos 26'i3)2 + sin^ 2i9i3 



Here, C is the same square root, which appears in the two neutrino matter formulas. 
Calculating the eigenvectors of M in order 0{a) gives: 



sin 2013 



a AsiT? 012 sirP 2613 



^/2C{A+C-cos2el3) 2Cy^2C^{-A+C+cos2ei3) 
a (1+A-C')sm2ei2siiiei3 

{l+A+c)^/2C{I+cP^^^^I2e^ 

sin 2013 a A sin^ 6>i2 sin^ 2gi3 



V2 



\ V2C(^XhC+^^^2e^ 2C^/2C^{A+C-cos2el3) I 

\ 



g cos 012 sin ^12 

A cos 013 

1 



a cos 012 sin 0i2 sin 0i3 



Acos2 013 
sin 2013 



+ 



a A sin 012 sin 2013 \ 



VaCC-A+C+cos 2013 )_ 26^26'^ (i+C-cos 2013) 
a (l+A-C) sin 20i2 sin 0i3 
(1+A+C')a/2C'(-A+C'+cos 2013) 

sin 2013 g A svc? 0i2 sin^ 20i3 



+ 0(g') 



\ V2C*(i+C'-cos20i3) 2(5^/2(52 (-i+C+cos 20i3) / 



(17a) 

(17b) 

(17c) 

(18a) 
(18b) 
(18c) 

(19) 



(20a) 



(20b) 



(20c) 



There is one major problem concerning the calculation of the eigenvalues and eigenvectors, which 
has to be addressed. Throughout the above series expansion A was assumed to be different from 
zero. This is important as the results given above do not hold for A = in which case a different 
series expansion in a would be obtained. This is a general and important fact. In principle, 
it is also possible to give results for small values of \A\, which, however, would fail for larger 
1^1 . The reason for this is that there are two different resonances occurring. One for A = a 
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(solar resonance) and one for A = cos 2^13 (atmospheric resonance). Each resonance produces a 
level-crossing of the eigenvalues. To describe both level-crossings, the correct expression for the 
eigenvalues are necessary. Being interested in approximative solutions, one has to distinguish the 
two above mentioned cases. In this work the focus is on the case \A\ > a, which is appropriate for 
neutrino beams above 1 GeV in matter densities of 2.8 g/cm'^ (Earth mantle) or more. However, 
one must not expect that the expressions for the mixing parameters in matter will show the 
correct convergence for A ^ 0. For Am'^i = 10~^ eV^ and 2.8 g/cm^ we find that A > a is valid 
for El, > 0.5 GeV. This lower bound on the neutrino energy decreases linearly with Am|^. 

That the results for the eigenvalues and eigenvectors obtained from the series expansion are not 
good at the resonance A ~ 1 is another point to mention. However, this does not have a crucial 
implication on the obtained results for the parameter mapping and oscillation probabilities. This 
issue will be discussed later, at the appropriate places. 



Construction of t/ 

It is now possible to construct U from the eigenvectors vi, V2, and V3. For this it is necessary 
to correctly identify the order and the signs of the eigenvectors. In order to avoid divergences 
in the expressions for the mixing angles, it is appropriate to change the order at the resonance 
A = cos 2^130 : 

jj ^ i {viV2V3)'^ for i< COS 26*13 ^^21) 

I {V3V2V1) for A > cos 2^13 

The second point is to bring U' to a form which is consistent with the standard parameterization. 
This is not trivial and has to be carried out carefully for each of the different cases. As an example, 
the case A < will be considered in detail: 

As the vacuum angle 623 was factored out from the beginning (eq. (^)), the matter induced change 
of this mixing angle 623 will be of order a. This can be also seen by looking at the (//, 3)-element 
of U. Furthermore, by looking at the (e, 2)-element, one finds that also 9i2 must be of order a. 
Considering this with the replacements §12 = as^^\ §23 = cts^^\ and S13 = sf'rj +as^^\ one obtains 
the following structure for U: 



/ . -(O)^(a) 

/ Ci3 OlC^sy Si3 



U 



-a(s 



(")+sS?4?) 1 ac^s^ \+Oia^). (22) 

(a) JO) , Ja) , 



\ -Si3 -a[Si2 Sl3 + S23 ) Ci3 

Then, sin^is and sin023 can be read off directly from C/e3, C/^3 and Ut3- 

sin 26*13 a A sin^ 6*12 sin^ 26* 



sin 013 = + -^^^-^ -13 ^ (23) 



2C{-A + C + cos 2^13) 2C^J2C'^{A + C - cos 2^13) 
. ^ (1+i-C) sin 2^12 sin gi3 ^ 

sm023 = a h C(a ) . (24) 



2(1 - ^ + C) cos 



^3 



* Another strategy would be to chose the order in such a way that in the limit \A\ —^0, the correct mixing matrix 
in vacuum is obtained. However, since the expressions for the eigenvectors and eigenvalues are not good in this limit, 
this is not a feasible solution here. 
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To find sin ^12, it is now useful to split off ^23- The rest tj = C/^(^23) U' should now be brought 
to the form 



.(0) .(a) 
CIS OiC\^s\2' Si3 

-q4? 1 



V 



'13 



.(a) ,{0) . 
-aSi2 Si3 Ci3 



The mixing angle 9i2 can then be read off from U'^i- 

a C sin 2^12 



smyi2 



(25) 



(26) 



A cos 6Ii3a/ 2C(-^ + C + cos 29 



13 J 



Parameter mapping 

Considering the correct ordering of the eigenvectors (eq. (|2T|)) and following the above described 
steps, one can determine the complete parameter mapping for all regions of the A parameter space. 
Comprising, one obtains the following expressions for the mixing parameters in matter: 

. ai sin 26*13 , a i sin^ 6*12 sin^ 26*13 
sintyi3 = — + ■j^^^^^^=^^^^^^= (z7aj 



2C{^A + C ± cos 2013) 2C'^^J2C{±A + C ^ cos 2^13) 
• Qi C sin 26*12 

sin0;^2 = 1^^^^^^=^^^^^^= (27b) 



1^1 cos 6*13 y 2C(=fA + C ± cos 26*13) 

. ni ■ , i sin 26*12 sin 6*13 COS 6*23 ^ 

smt/oq = smcy23 + «coso (27c 

±l + C7^Acos2% 

, / cos (5 2I sin 26*12 sin 6*13 \ 

smo = smo 1 — a (27d) 

\^ tan 26*23 ±1 + C =F A cos 26*13 / 

Here, in the expressions with choices for the sign, the upper sign holds for A < cos 2^13 and the 
lower sign holds for A > cos 2^13. Higher orders than 0{a) are omitted. To take into account 
also ^23 and 5, which were factored out at the beginning, the equations (p^a-d) were applied. The 
expansion of sin (5' given here does not hold for ^23 ~^ 0- 

From this parameter mapping it is possible to derive the following quantities: 

. ^ sin2 2% 2i(-i + cos2%)sin2 0i2sin2 2% 
sm 26*10 = ha (28a) 

2(7 sin 26*12 . ,^ 

sin 20^2 = a \ (28b) 

|i| cos 0i3y 2(7(^i + C ± cos 2^13) 

■ na' ■ ofl , A 2^ ^'^^ ^^12 sin 6*13 cos 26*23 ^ 

sm2t^no = sm 26^23 + a cos (28c) 

±l + C7^^cos2% 



For the mass squared differences one obtains: 



2'a 2'a 2'\ f( Amo, Attio, Amf) for ^< cos 26*13 /<ir.\ 
^ ' \ (-Amf,-Ami,-Am§) for yl> cos 26*13 
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with 



Amf := A(A3 - A; 



2 



-(1 + A + C)A - aA I COS 612 ^ I , (30a) 



Amg := A(A2 - Ai) 



1/ J /vx. . 2n (1 + (7-1008 2(913)31112(912 1 , 

= -(-l-^ + C)A + aA ( cos20i2-^ 2Cr^ ^J' 

Am^' := A(A3-Ai) 

= (7A + a^(-^ + -^-°!^^^^)-^"^^^^ (30c) 

(y 

Looking at the expressions for the mixing angles in matter, one obtains the following interesting 
statements: 

sin^ 26>;3 

In leading order, one finds the well-known resonant behavior of ^13 familiar from two neutrino 
oscillation as MSW-resonance. The order a correction to this leading result is suppressed by two 
powers of ^13, and hence, is negligible small. A careful study of the correction indeed shows that 
it is small and only important if precise results are to be obtained. The expressions for ^13 do 
not show divergences for |^| — > and the vacuum limit is correctly described. Comparison with 
numerical results shows an excellent agreement even for \A\ < a. 

sin 20^3 

In leading order, the mixing angle ^23 is equal to the vacuum mixing angle sin 2023- The order a 
correction is double suppressed by ^13 and by cos 2^23 (when ^23 is close to 7r/4). Its proportionality 
to cos 6 is caused by the mixing of the CP-phase 6 with the 0{a) correction of ^23 (^q- (|T5|c)). The 
expression for ^23 shows the correct behavior for \A\ and numerical results are consistent also 
for \A\ < a. 



sin 20^2 

The quantity sin 20^2 is of order a. For a ^ it does not reproduce the vacuum parameter ^12 • 
But this is not difficult to understand. For = 0, the first term in the Hamiltonian (eq. (^)) is 
invariant under rotations in the 12-subspace. This reflects the fact that for a = the solar mixing 
angle does not influence the oscillation probabilities and could in principle be chosen arbitrarily. 
Interesting here is that sin 20^2) even for large values of \A\, is proportional to a. In leading order 
of ^13 one finds that sin 20^2 = sin 2^i2/|>l|. There appears a divergence for \A\ 0. The 
result is unphysical for \A\ < a, which reflects the problem that the level crossing at the solar 
resonance is not correctly described. Since |^| is proportional to the neutrino energy E,y, sin2^^2 
is suppressed not only by the mass hierarchy, but also by large neutrino energies. 
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CP-phase S 



The correction to the CP-phase 5 in matter is triple suppressed by the mass hierarchy a, 613, 
and tan~^2023- For sin^ 2^23 = Ij the CP-phase 6 is not changed (in order a). The invariance 
of sin 5 sin 2^23 under variations of the matter density p (eq. (^d)) is an exact result, which is 
independent from the approximations made. 



4 CP-violation: Jcp in matter 



rim 

violating effects. In vacuum, it is given by 



From the vacuum case it is known that the quantity Jqp = Ini J-J^ drives the strength of CP- 



8 Jcp = sin 5 cos ^13 sin 2^12 sin 2013 sin 26*23 • (31) 
Application of the parameter mapping (eqs. (|27|)) gives J^p in matter: 
sin 5' cos 6*13 sin 26*12 sin 26*13 sin 26*23 = 

sin 6 cos 013 sin 20i2sin20i3sin2023 + O(a^) . (32) 



|A|Ccos2 0i3 
One thus finds the important and simple result 

•^cp = TTTs; TT— -^cp • (33) 

|A|Ccos2 0i3 

Applying this result to 

JcpAm'i2Am3iAm32 = JcpA^a + 0(0^) , (34) 

the Harrison- Scott invariance jQpAmi2Am3iAm32 = JcpAmi2Am3iAm32 § can be verified. 

It is important to notice that also in matter all CP-violating effects are proportional to the mass 
hierarchy a. In vacuum, the suppression of CP-effects through the mass hierarchy is obtained from 
the smallness of the solar mass splitting, which is aA. In matter, the mass hierarchy is lifted, but 
the mass hierarchy suppression is retrieved in sin20'i2, which is proportional to a, and thus, leads 
to a mass hierarchy suppression of J^p- 

Another interesting point to notice is the factor 1/(7, which leads to an MSW-like resonant en- 
hancement of Jqp in matter. It can thus be expected that the CP-terms PsinS and PcosS will 
benefit from the MSW-resonance in the same way as the leading two neutrino term Pq does. 



5 The Ue — > appearance probability 

Having presented the parameter mapping in matter, it is now possible to start from the ordinary 
vacuum expressions (eq. (|l|)) in order to derive the oscillation probabilities in matter. The Jfp as 
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series expansion in a take the following shape: 



Im Jll^' 



Im Ji'^g^' 



Re Jit' 



ReJ,t' 



Im Jl^' 



cos 5' sin cos^ ^13 sin cos ( 



sm 



9' 



(35a) 



-sin^ cos^ ^23 + 
— cos b' sin 6*^2 cos^ Q'y>^ sin 6*^3 cos ^23 sin ^23 
-sin2 20'i3sin2 023 + C(a^) (35b) 
cos b' sin 0^2 cos^ ^^3 sin ^^3 cos ^23 sin ^23 ^0(0?) (35c) 
cos b' sin 0^2 cos^ 6*13 sin ^^3 cos ^23 sin 6*23 ^0{o?) (35d) 



Even though in general the calculations were performed only up to order a, a closer look at o?- 
terms proves to be important. Each second term of Re jf in eq. (|35| ) is of order . Since Q'y2 
is not suppressed by ^13, these terms give a non- negligible contribution to the overall oscillation 
probability. This order c? sin*^ ^13 contribution, which will be identified with the P3-term in vacuum 
(eqs. dll)) is important for small values of Q\z- It is possible to show without explicit calculation 
of all order a^-terms of the parameter mapping that no further terms of this kind exist. All other 
Qf^-terms in the oscillation probability will at least be suppressed by one power of ^13. 

Inserting the expression for the mixing parameters in matter together with the abbreviation A = 
gives the following list of terms contributing to the oscillation probability P(i^e — *■ v^: 



sin^ 29 



sm 



723" 



13 



fj2 



sin2(AC') 



R 



sin 5 



1 sin b cos ^13 sin 2^12 sin 2^13 sin 29 



-a 



23 



AC cos 



1 cos b COS 6*13 sin 26*12 sin 26*13 sin 2^23 



AC cos 6*^3 



(36a) 



sin((7A) cos(C'A) - cos((l + i)A) (36b) 



sin(CA) sin((l +A)A) ^sin((:7A) (36c) 



1 - A cos 26* 



-a 



13 



a 



C3 

2i(-i + cos 26*13) 
^Fl + Ciicos 26*13 



a 



2CMcos2 ^13 

2(7 cos^ 6*23 sin^ 26*12 



sin^ 6*12 sin^ 26*13 sin^ 6*23A sin(2AC') 
sin^ 6*12 sin^ 26*13 sin^ 6*23 sin2(AC') 
cos 013 sin 2012 sin 2^13 sin 2023 sin^ (AC) 



(36d) 
(36e) 
(36f) 

by flipping the sign 

of the -Psin5 term. In all expressions with two possibilities for the sign, the upper sign is valid for 
A < cos 2013 and the lower sign is valid for A > cos20i3. The A-dependent pre- factors of Pi,P2, 



A^ cos^ 013 (=p A + C lb cos 20 



13 j 



sin^ -(l + iTC)A 



The probability P{ve 



can be obtained from the probability P{ve 
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and P3 expanded in ^13 give: 



C3 {A - 1)2 

2i(-i + cos26li3) _ 2A 

^4 ~ ~(i-l)3 
Tl + C±icos2013 2 ^ 



2C72Acos2 % 
2(7 



1+O(0?3) 



cos2 6li3(=F^ + C±cos 2%) 

Thus, Pi is quadratic in sin0i3 and P2 even of third order in ^13. Therefore, Pi and P2 are 
neghgibly small compared to Psin5 and PqosS- The term P3 is important, since it is the only term, 
which is not suppressed by ^13. It was stated before that in some cases the expressions for the 
eigenvalues and eigenvectors are not good at the resonance A = cos 2^13. This problem stems from 
the second order in ^13. On the level of probabilities, this deficiency is small and only visible in 
the Pcos(5-term for large values of ^13. It turns out that neglecting the subleading terms, which are 
the source of this problem, gives very accurate results also for A = cos 2^13. This modification can 
be applied to both the Pcos5-term and the Psin5-term: 

sin (5 cos % sin 26^12 sin 26*13 sin 2023 ■ K ■ a ■ a a ro^ \ 

Psm5 = a — :^ smGAsmAsmAA , (37a) 

AC cos 6*^3 

cos (5 cos 6*13 sin 26*12 sin 26*13 sin 26*23 . a; A • ;i a /Q^7U^ 

Pcos5 = ce — sm G A cos A sm . (.-J'b) 

AC cos 6*f3 

Neglecting all subleading terms in ^13, the relevant terms Pq, Psin(5) PcosS^ and P3 take the following 
simple shapes: 

Po = sin2 023S^sin2((i-l)A) (38a) 

P^^s = a ^ sin(A) sin(iA) sin((l - A) A) (38b) 

^(1 — ^4) 

cos <5 cos 6*13 sin 26*12 sin 26*13 sin 26*23 / a\ ■ f a a\ ■ f ft a\a\ fia \ 

PcosfS = ct x cos(A) sm(AA) sm((l — ^)A) (38cj 

^4(1 — A) 

2 Cos^ 6*23 sin^ 26*12 . 2/Ja^ /QQ^^ 

P3 = a J- sm [AA) (38d) 

^2 

It is evident that in the limit of small baselines, A — > 0, these expressions converge to the results in 
vacuum (eqs. (|3|a-d)). A numerical study shows that the precision loss of eqs. (p^a-d) compared to 
eqs. (^6^-f) is only relevant for the largest allowed values of sin2 2^13 near the CHOOZ-bound (0.1). 
The precision loss is mainly caused by the approximations made in Pq . The term P3 contributes to 
the overall probability only for small ^13, and hence, does not suffer an appreciable accuracy-loss 
in the form given in eq. (|38|d). Figure |l] shows a comparison of the analytic results obtained here 
with the results obtained from a numerical study. Note that the combined contributions from 
eq. (|3^a), eqs. (^a,b) and eq. (|38|d) are identical to the result obtained by Cervera et al. [^] (eq. 
(16)). A similar approach has been discussed in ref. |jl^. However, eq. (16) therein does not cover 
the case of very small ^13, since it does not include order {Arn^i/ Am'^^)'^ corrections. 
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Figure 1: Analytical results (dashed and dotted lines) compared to numerical results (solid line) for the 
oscillation probability P(i^e ^ i^^) in matter (2.8 g/cm'^) as function of the neutrino energy. Negative 
energies correspond to anti-neutrinos. The dashed line uses the expressions ^6|a,b,c,f. The dotted line 
was obtained from equations |3^a,b,c,d. The calculation was performed for the baseline L=7000 km with 
3 = 0, bimaximal mixing and three values of sin^ 26*13 (0-1, 0.01, 0.001). The squared mass differences are 

^31 



Am?i = 3.2 • 10^3 and Am^, = 1 • lO^^eV^ 



6 Applications 

6.1 Validity region of the low L/E,y approximation in matter 

Frequently, the low L/Eu limit is used to simplify complex calculations or derive power laws for 
neutrino rates. In vacuum, it is well-known that this approximation is valid for 

A<1 ^ E,> 4.0GeV ( ^""^ , ] ( ^] . (39) 

V 3.2- 10-3 eVV V 1000 kmy ^ ^ 



With the use of eqs. (pqa-d) , it is possible to extend this argument to the presence of matter. Note 
that in the oscillatory terms, which are linearized in the small A approximation, there now also 
appear the terms ^A, which must be small. In this product, the dependences on the energy E^^ 
and the mass squared difference Am|| cancel. Hence, in addition to relation (^), a direct limit 
on the baseline L, which only depends on the matter density p is obtained: 



AA<1 => L< 3700 km ^ 3 . (40) 

\2.8 g/cm'^y 

6.2 CP-asymmetry in matter at small L/E,^ 

CP-violation studies frequently focus on the fundamental quantity called CP-asymmetry Aqp: 
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In vacuum, being proportional to sin (5, ^cp is a direct measure for intrinsic CP-violation. Since 
^CP is a ratio of probabilities, it has the important advantage that, on the level of rates, systematic 
experimental uncertainties to a large degree cancel out. However, matter effects also create fake 
CP-asymmetry, which spoils measurements of the intrinsic CP-violation induced by 5. The problem 
to distinguish these two different sources of CP-violation is often called the "disentanglement 
problem". In a typical long baseline neutrino experiment, the strength of matter induced CP- 
effects reaches the strength of intrinsic CP-effects at baselines around 1000 km. 

Using the above derived approximative solutions for the appearance probability P(fe 
possible to calculate the small A limit of vlcp- For bimaximal mixing (^23 
given by 



it is 



5*12 = vr/4) ^cp is 



^CP 



2A sin 26*13 cos 6*13 {aAA cos 5 - SasinS + 2AA sin 6*13) 



1 



(42) 



3(a2 + 2a cos 6*13 sin 26*13 + sin^ 26*13) 

The approximation is valid in the regime given by eqs. (|3^) and (^0|). This limit is helpful to 
describe the behavior of Aqp for higher neutrino energies at not too long baselines. It is interesting 
to notice that in principle the leading contribution to ^cp in ^ has its origin in the sin 5 term. At 
first sight, this would suggest to distinguish this intrinsic contribution from matter contribution 
of order by the energy dependence of Acp. However, taking into account that A itself is 
proportional to E^, it turns out that all terms in eq. ( |4^ ) have the same energy dependence 
To summarize: In leading order in A, the CP-asymmetry in matter is proportional to l/Ey. The 



o 

" 2 

o 
< 



sin^ 26*1 3 =0.001 




^/2 nl2 
CP-^hase 6 



-fr/2 nil 
CP-phase 6 



-n -n/2 n/2 n 
CP-phase 6 



Figure 2: Dependence of the high energy limit of the CP-asymmetry on the CP-phase 5 for bimaximal 
mixing. On the ordinate is plotted the value of E^Aqp in GeV, which should be energy independent in the 
low L/Ei, approximation. The solar mass splitting was chosen at the upper edge of the LMA-MSW solution 
Am2i = 1 • 10^^ eV^ and the atmospheric mass splitting was varied in the Super-Kamiokande allowed 90% 
confidence interval 3.2 • 10""^ < Arn^^ < 3.6 • 10^"^ eV^. The calculation was performed for a baseline of 
1000 km. 



coefficient, which describes the l/£',/-energy dependence of Aqp for high energies is sensitive to 
both, matter effects A and intrinsic CP-effects from 5. At high energies, the quantity E^Aqp 
is predicted to be constant in the energy spectrum and this characteristical quantity could give 
direct access to the CP-phase 6. This is demonstrated in fig. |2|, which shows the value of E^Aqp 
as function of the CP-phase 6 at different values of sin^ 2^13. Since E^Aqp does not vary with the 
energy, this simple analysis is to a good approximation independent from the energy distribution 
of the neutrino beam. It is of course questionable if, in a real experiment, in the constant regime 
of E^Aqp, there are enough neutrino events to measure. Also this method cannot replace a full 
and detailed statistical analysis of the complete neutrino energy spectrum. 
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6.3 Strength of the CP-terms Psin<5 and PcosS 



The two subleading terms (|36|b) and (pqc) currently raise considerable interest as they contain 
information about the CP-phase 6 of the neutrino sector. Today, much effort is spent on the study 
of CP-violating effects in neutrino oscillation experiments p^ . One can try a simple approach to 
this problem by using the here obtained analytic results. It would, for example, be interesting to 
know, how strong the information on 6 inherent to the appearance oscillation probability is. To 
quantify this, one can look at the relative magnitude of \PsmS + -fcos^l compared to the statistical 
fluctuations \/Pq + in the background signal (provided the errors are Gaussian). To obtain 
statistical meaningful numbers, the estimation should be performed at the level of event rates 
expected in a real experiment, e.g. a neutrino factory long baseline experiment. Typically, flux 
times cross sections of a neutrino factory beam |12] scales like E'^/LP' . A neutrino factory of 20 GeV 



muon energy and lO^'^ useful muon decays per year produces 54800 f^-events in a 10 kton detector 
at 1000 km distance (assuming measurements in the appearance channel). As a statistical estimate 
the following ratio could be chosen: 




E in GeV E in GeV E in GeV 

Figure 3: Icr, 2(T, 3ct and 4(T contour lines of the quantity S (eq. (^)) in the L — E^^ parameter plane. 
Light shading indicates no signal and dark shading indicates strong signal. The left plot studies only the 
^siniS term. The plot in the middle displays the strength of the PcosS term. The right plot, which combines 
both terms should give the best approximation to more complex studies. Note that no energy spectrum was 
used in this crude model. The calculations were performed with (5 = 7r/2 (left), 5 — Q (middle), 5 — tt/A 
(right), bimaximal mixing, and sin^ 20i3 = 0.01. The mass squared differences are ^m^^ — 3.2 • lO^'^eV^ 
and Amii = 1 • IQ-^eV^ 



The value of S gives the number of standard deviations ("cr's") at which the CP-signal is distinct 
from the "background" . Figures |3| show the contour lines Icr, 2cr, Scr and 4a" of S in the L — 
parameter plane. The plots were produced with a running average matter density matched to the 
baseline L. It is interesting to note that in most of the L — Efj_ parameter space, there is no obvious 
decrease of the statistical sensitivity to CP-effects for increasing beam energy E"^ as often quoted 
in the literature. To study this point in more detail, it is helpful to derive the low L/Ey (eq. 
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scaling laws for 5 in the cases sin 5 = 1 and cos 5 = 1: 

SsinS ' j= and Scoss ~ • (44) 

Indeed, for the PsimS-term, the statistical sensitivity should decrease like Xj ^E^. However, the 
validity-region of the low LjEy approximation, according to eq. (|39|), is Ey > (4, 12,20) GeV for 
L = (1000, 3000, 5000) km. In the left plot of fig. ^ it can be seen that roughly at these energies, 
S shows a plateau where its maximal value is reached. The argument in favor of small energies 
thus only holds for very small baselines around 1000 km and smaller. The sensitivity to the Pcosb- 
term increases like y^E^j_. Hence, in the case of large cos (5, high beam energies are favored to 
extract information on the CP-phase 5. In conclusion, the difference of the result presented here 
and statements being found in the literature has two sources. First, usually only the explicitly 
CP-violating part PsinS of the oscillation probability is assumed to give the CP-signalQ. Second, 
the high energy approximation to the oscillation probabilities is often applied without careful 
consideration of its validity region. 



7 Conclusions 

The purpose of this work was to find approximate analytic expressions for the neutrino mixing pa- 
rameters and oscillation probabilities in the presence of matter. It was stated that being interested 
in approximate solutions it is difficult to describe both the solar and the atmospheric resonance at 
the same time. Therefore, this work is restricted to energies above the solar resonance according 
to: 

|i|>H^^„>0.«Gev(^)(^). ,45) 

For this regime, the complete parameter mapping (eqs. (p7|)) was given as series expansion in the 
small mass hierarchy parameter a = Am2^/Am|^. It was shown, that the change of the CP-phase 
5 in matter is triple suppressed by the mass hierarchy, the mixing angle ^13 and by ^23 being close 
to maximal. Furthermore, it was shown that in order Amg^/Aml^^, the relevant contribution to 
the parameter mapping is the correction of 612 in matter. The derived parameter mapping was 
used to compute the P{ve ^^i) appearance oscillation probability in matter. Effort was made to 
find simple solutions, which hold over a wide parameter range and are easy to compare with the 
results known from vacuum oscillation. An answer, which in the author's point of view fulfills all 

^ Frequently, the need for explicit detection of an asymmetry between the two CP-conjugated channels is stressed 
and matter effects are considered as background, which prevents such measurements. The attitude taken here is, 
however, different: The goal of any experiment is the limitation of the allowed parameter space for 5, which does 
not necessarily presume the detection of explicit CP-violation. Hence, the Pcos 5 contribution has the same status as 
the Psin(5-term and matter effects have to be included in the theoretical model, which is fitted to the experimental 
data. 
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these requirements is the fohowing set of terms (eqs. (|38|)) contributing to P(z/e — > z/^): 
= sin2 023^5^ sin2((i - 1)A) , 

sin COS 6*13 sin 26*12 sin 2013 sin 26^23 ■ r k\ ■ / a K\ ■ t ^^ 

PsinS = a J- J- sm(A) sm(^A) sm((l - ^)A) , 

^(1 - A) 

cos (5 cos 6*13 sin 26*12 sin 2013 sin 26*23 r a \ ■ ^ a a\ ■ r A\/i\ 

PcosS = Oi cos(A) sm(ylA) sm((l - ^)A) , 

^(1-^) 

2 cos^ ^23 sin^ 2012 . 2/^A^ 
Pa = J- sm^(ylA) 

A^ 

with A = /\ml^L / {AE^) and A = Aj I^m\^ = 2VE^/Aml^. This gives quaUtatively good results 
for basehnes at which the oscihation over the small (solar) mass squared difference can safely be 
linearized^: 

oA<l^L<8000 k,n(^)(ll^). (46, 

To obtain high precision results for large values of 0i3 , it is recommended not to neglect subleading 
013 effects. The corresponding terms to P{i'e — > i^fj.) are given by eqs. (36a,b,c,f). Results for the 
anti- neutrino channel are always obtained by flipping the signs of Psin<5 and A. 

Using the derived approximations to the oscillation probability, it was shown that from relation 



(40) a stringent limit on the baseline L can be derived, up to which the small L/E,y approximation 
in matter is valid. Then, using this approximation, an expression for the CP-asymmetry Aqp in 
matter was given, which demonstrates that, for high neutrino energies, Aqp is decreasing pro- 
portional to 1/E^. It was proposed that measuring this energy-dependence could help to obtain 
information on the CP-phase 6. Last, it was demonstrated that estimations on the experimental 
sensitivity to the CP-terms in P^i'e i^^) can be given. The here obtained results do not favor 
low neutrino energies for the CP-violation search. The reason for the discrepancy between this 
result and statements, which can presently be found in the literature, were discussed. These topics 
were discussed only briefly and mainly serve as demonstrations of the applicability of the derived 
formulas. 
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Of course, it is also possible to give results, which are not limited by this baseline restriction. However, this 
approximation is very helpful to obtain simple results. 
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